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ABSTRACT 

We present self-consistent global, steady-state MHD models and synthetic optically thin synchrotron 
emission maps for the jet of M87. The model consist of two distinct zones: an inner relativistic outflow, 
which we identify with the observed jet, and an outer cold disk-wind. While the former does not self- 
collimate efficiently due to its high effective inertia, the latter fulfills all the conditions for efficient 
collimation by the magneto-centrifugal mechanism. Given the right balance between the effective 
inertia of the inner flow and the collimation efficiency of the outer disk wind, the relativistic flow 
is magnetically confined into a well collimated beam and matches the measurements of the opening 
angle of M87 over several orders of magnitude in spatial extent. The synthetic synchrotron maps 
reproduce the morphological structure of the jet of M87, i.e. center-bright profiles near the core and 
limb-bright profiles away from the core. At the same time, they also show a local increase of brightness 
at some distance along the axis associated to a recollimation shock in the MHD model. Its location 
coincides with the position of the optical knot HST-1. In addition our best fitting model is consistent 
with a number of observational constraints such as the magnetic field in the knot HST-1, and the 
jet-to-counterjet brightness ratio. 

Subject headings: MHD - methods: numerical - radio continuum: galaxies - galaxies: jets - galaxies: 
individual: M87 



1. INTRODUCTION 

Since its discovery by iCurtid (|1918t ) the jet of M87 is 
the classical prototyp e for extragalactic jets. Due to its 
proximity at 16 Mpc (IWhitmore et alj|1995t iMacri et ail 



/I PC 

119991 : iTonrv et al.N2001h M87 is one of the closest radio 
galaxies, which allows present VLBI instruments to re- 
solve the transversal structure of the jet. Therefore, it 
is an ideal candidate for testing specific jet formation 
mechanisms. The jet and its hot spots have been sys- 
tematically studied across the electromagnetic spectrum 
from the radio to X-rays, both, with ground-base d obser- 
vatio ns and from satellites (for a review see e.g. iBirettal 
1996). The initial opening angle is approximately 60° on 
scales of about 0.04 pc and decreases rapid ly until reach- 
ing 1 0° at a distance of 4 pc from the core (Biretta et al.l 
2002). These observations suggest that the jet of M87 is 
rather slowly collimated across a length of several parsec. 

The prevailing paradigm for jet formation and 
collimation is magnetic self-collimation by the 

Electronic address: jgracia@cp.dias.ie 



iBlandford k Pavnd {1982) mechanism. This view 
is supported by observations, which are consistent with 
the presence of a non-van i shing t oroid al magnetic field 
component (|Asada et aTl 120021 120081: iGabuzda et al.l 



12004 IZavala k Taylor! 120051 : iGomez et al.ll2008D . How- 
ever, relativistic effects have been shown to decrease 
the collimation efficiency, i.e. for a given magnetic 
field configuration at the base of the jet (or rotator 
efficiency) , the fraction of total mass- and magnetic flux 
that is asymptotically cylindrically collimated is lower 
for a relativistic flow than fo r a non-relativistic flow 
(iBogovalov k Tsinganosl 119991 : iTsinganos k Bogovalovl 
2000]). Not only is the mass flux fraction lower, but the 
final opening angle is larger for relativistic outflows due 
to the decollimating effect of t he electric field and the 
effective inertia of the plasma (Bogovalovl |2001[ ) . which 
both counter-act the pinching by the toroidal magnetic 
field. However, this is strictly true only for initially 
radial magnetic field structures, as opposed to extended 
magnetic field co n figurat ions like disk- winds. S ee, e.g . 
iFendt k Memolal (j200lh : iVlahakis k Koffigl |2004); 
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iKomissarov et all ([2007) for efficiently collimating MHD 
disk- wind models. 

of papers 
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series ot pa pers Itiracia, et alj 
iTsinganos fc Bogovalovl (|2005l . [2002) suggested a 
steady-state two-component MHD model. The model 
consists of an inner relativistic outflow, which is identi- 
fied with the observed jet, and an outer non-relativistic 
disk- wind. While the inner relativistic jet is not expected 
to collimate well through magnetic self-collimation, the 
very same process operates efficiently in the outer 
non-relativistic disk-wind. It is expected, that at least 
for a part of the available parameter space, collimation 
in the outer disk- wind is so efficient, that it might 
resist the decollimating inertia of the inner relativistic 
plasma and channel t he jet into a narrowly confined 
beam. iGracia et al.l (|2005f ) have shown, that such 
two-component models could easily account for the 
narrow appearance of the beam of extragalactic jets 
by reproducing the observational measurement of the 
opening angle distributi on as a func t ion o f angular 
distance from the core bv lBiretta et~aTI (|2002fh bmce m 
these models the jet, i.e. the relativistic inner outflow, 
is strictly speaking not magnetically self-collimated, 
but rather confined by the outer disk-wind, the authors 
prefer to ta lk of collimation by magnetic confinement. 

However, IGracia et al.1 (|2005l ) could not fit the open- 
ing angle distribution with a unique MHD model. In- 
stead various sets of parameters reproduce the observa- 
tions with similar accuracy. 

In deriving the opening angle of their model, 
IGracia et alj (|2005f ) made a simplifying but crucial as- 
sumption. They identified the observed jet with the in- 
ner relativistic outflow of their two-component model. 
More specifically, the boundary of the jet was assumed 
to coincide with the shape of a specific fieldline the 
one fieldline, which separates the relativistic inner region 
from the non-relativistic outer disk-wind at the base, or 
the launching surface, of the outflow. So, the observed 
opening angle was strictly speaking fitted by the shape 
of a single magnetic fieldline. 

However, observations do not measure the plasma state 
directly, i.e. in terms of velocity, temperature or mag- 
netic field strength. Instead, they register photon flux as 
a function of position on the plane of the sky. As such, 
from an observational point of view, the width of the jet 
is defined by the emission dropping below the detection 
limit or a small fraction of the luminosity of the ridge line 
of the jet. So the question is how well did IGracia et al.1 
(2005) measure the width of the jet in terms of observa- 
tional quantities? Or more generally - Do MHD models 
explain the appearance of AGN jets? The purpose of this 
paper is to answer exactly this question by adopting the 
point of view of an observer. Assuming, that the main ra- 
diation mechanism is synchrotron emission, we translate 
the steady-state MHD model into a synthetic emission 
map and measure the width of the jet using only these 
data. 

The outline of the paper is as follows. In the following 
section [5] we summarize the two-component MHD model 
for extragalactic jets and discuss some of its properties 
relevant to this work. In section [3] we present compact 
expressions for the calculation of synchrotron emission 
and apply these to our numerical MHD models for the 
jet of M87. Finally, we discuss our results and draw some 
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Fig. 1. — Illustrative sketch of the model (see text). 
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Fig. 2. — Comparison of the opening angle calculated from 
our MHD models and the observational data for M87. The lines 
show opening angle profiles as given by the separating fieldline \P Q 
for model A (solid line), i.e. the best model, and model D (dot- 
dashed), which best reproduces the data based on \l/ a . The models 
B ((dashed)) and C (dotted) are shown for completeness. The data 
points marked by filled circles were taken into account in the fit- 
ting procedure. The innermost and outermost measurement (open 
circles) were disregarded. 

conclusions. 

2. MHD MODEL 
We adopt the model and notation of IGracia et al.l 



(2005) and refer the reader to that paper for details. A 
simple illustration of the model is shown in Figure [TJ It 
consist of two distinct zones; an inner outflow, which 
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TABLE 1 

Parameters and description of models. 



name 


a[»] 


Tj [mc 2 


r .i 


Wj [r g /c] 


r d 


Bo 


remarks 


model A 


16.15 


2.85 


2.74 


3.0 


1.02 


1.1 


best overall model 


model B 


14.45 


3.0 


2.74 


3.0 


1.02 


0.9 


best fit to opening angle based on synchrotron map 


model C 


24.0 


3.0 


1.8 


2.7 


1.02 


1.0 


highest limb-brightening and brightest knot HST-1 


model D 


14.4 


3.45 


3.26 


3.0 


1.01 


1.27 


best fit to opening angle based on field line 


The mass 


flux rate of 10 24 


gs 1 is a parameter common to all models. 
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Fig. 3. — Lorentz factor along the flow for model A. Contour lev- 
els are shown at T = 2,4, .., 12. The relativistic inner jet initially 
has Lorentz factors r ~ 2 — 3, but accelerates up to T ~ 10 — 12. 
Strong gradients may be seen across the jet and near the recolli- 
mation shock, if present. 

is dominated by relativistic dynamics, and an outer 
non-relativistic outflow. Both outflows originate from a 
spherical launching surface located at a distance rrj from 
the black hole. The launching surface is threaded by a 
helical magnetic field; the poloidal component is initially 
perpendicular to the launching surface. The two zones 
are separated by a specific fieldline ^f a , where 2a is the 
initial angular width, or opening angle, of the inner rel- 
ativistic outflow. In the following, we will refer to these 
two distinct zones by relativistic jet and (non-relativistic) 
disk-wind, respectively. 

This two-component model is motivated by a similar 
two-component structure of the und erlying accretion flow 
consi sting of an outer standard disk (jShakura fc Sunvaevl 
fl97l and an inner hot plasma, which c ould be either 
an advection dominated accretion flow (|Naravan & Yi 
[1991 iPeitz fc App]|[l997t iGracia et aLll200.3h . or the fi- 
nal plunging region near the black hole, where relativis- 
tic dynamics dominates through, e.g., frame-dragging, or 
the Blandford-Znajek process. 

We impose two different sets of boundary conditions 
in the two distinct zones along the launching surface. If 



the launching surface is located beyond the fast magne- 
tosonic surface, i.e. in the hyperbolic MHD regime, the 
steady-state problem reduces to an initial value Cauchy- 
type problem and the steady-state equations can be in- 
tegrated directly in terms of conserved integrals o f mo- 
tion as described in iTsinganos fc Bogovalovl(|2002l ). We 
stress, that beyond the launching surface we solve the 
axisymmetric steady-state problem self-consistently, in- 
cluding the magnetic field structure, as a function of the 
boundary values alone. However, we do not solve the 
problem inside the launching surface, which is a much 
more complicated exercise. A self-consistent solution of 
the full problem needs to take into account the dynam- 
ics of the accretion flow, something which is beyond the 
scope of this paper. 

This procedure yields a set of quantities as a function 
of space in the comoving frame of the jet. The sepa- 
rating fieldline ty a perfectly divides the relativistic from 
the non-relativistic outflow. Inside of the plasma is 
highly relativistic, both in terms of its bulk Lorentz fac- 
tor, r > 1, and of its thermal energy, T > mc 2 , while 
outside the plasma is cold, T <C mc 2 , and moving at 
non-relativistic speeds, T = 1. Also, the magnetic field 
strength peaks close to the separating fieldline, where 
fieldlines (and poloidal flowlines) of the inner outflow are 
strongly compressed laterally and confined to a narrow 
sheet by the fieldlines of the outer disk- wind thus forming 
a natural interface be tween both zones. 

IGracia et al.1 (|2005l ) used the shape of the separat- 
ing fieldline t o fit the observed opening angle 
(jBiretta et al l [20021 as a function of distance from the 
core, however, without taking into account projection ef- 
fects, i.e. they assumed implicitly, that the jet of M87 
was in the plane of the sky. There is an on-going dis- 
cussion on the inclination angle of t he M87 jet (see e.g . 
Owen et all 119891 : iReid et all Il98fl iBiretta et all 11999b 



Lv et al.ll2007f L In this paper we assume, that the jet of 
M87 is oriented at an angle 8\ os = 40° from the line-of- 
sight as a compromise of values discussed in the litera- 
ture. 

Unfortunately, the parameter space of the MHD model 
is degenerated in the sense, that very different sets of pa- 
rameters yield equally good fits or even almost identical 
opening angle distributions. Then, in order to further 
constrain the model parameters we shall invoke the radi- 
ation signatures of each model and compare them with 
the corresponding observations. In this way, we may pin 
down a small number of acceptable models which simul- 
taneously satisfy the constraints of the MHD model and 
also reproduce the observed distribution of the emitted 
radiation. 

We have run 2600 axisymmetric steady-state MHD 
models and calculated synthetic synchrotron maps for 
them. Here we discuss the four models that evaluate 
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Fig. 4. — Comparison of the normalized synchrotron intensity 
along the jet axis calculated from model A (solid lines) and model 
C (dashed lines), respectively. The dashed lines indicate power- 
laws with index -1.5 and -2.5, respectively. 

best under different criteria. See table [1] for their param- 
eters. Figure [2] compares the opening angle as defined by 
the separating fieldline projected under an angle of 
#ios = 40°, with the observational measurements. Note 
that all our models fail to fit the jet width at large dis- 
tances close to the optical knot A at ~ 900 pc and also 
at small distances close to the jet's origin, where there 
are some uncertainties, as is discussed in the last section 
of the paper. Models with large initial opening angle 
matching the first data point as e.g. model C, have dif- 
ficulties reproducing the opening angle distribution. It 
is in general very difficult to make the curve more con- 
cave. Models with initial opening angles falling slightly 
short of the measured value (as models A, B, and D) 
may easily reproduce the rest of the measurement and 
yield quite good fits. The quality of agreement between 
the observed opening angle distribution and the theo- 
retical curve is measured by a simple x 2 ■• i- e - the sum 
of squares of the difference over the available observa- 
tional data points. Model D best fits the opening angle 
at x 2 = 2.9, but does not reproduce the morphological 
structure of the radiomaps particularly well, as will be 
discussed in the next section. However, the overall best 
model A reproduces the morphological structure and still 
fits the opening angle data reasonably well with x 2 = 7.8. 
Note, that both models show a clear kink at ~ 100 pc due 
to recollimation towards the axis. 

Typically, models that fit the observed opening angle 
well are moderately relativistic, both in terms of the ini- 
tial outflow velocity Tj ~ 2 — 3 and the initial plasma in- 
ternal energy 2] ~ 3mc 2 . In these models the plasma ve- 
locity increases along the flow to values up to T ~ 5 — 10. 
However, the plasma may decelerate and reaccelerate 
sharply at the recollimation shock. Strong gradients of 
the plasma velocity may generally be present across the 
jet, as seen in Figure. [3] It is therefore very difficult to 
assign a typical Lorentz factor to the whole jet. 

3. SYNCHROTRON MAPS 



normalized intensity 
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Fig. 5. — Convolved synthetic synchrotron map. To increase 
contrast, the map has been divided by the trend along the jet axis 
(Figure [4}. The two lines near the edge of the jet indicate the jet 
width as defined by the separating fieldline (inner black line) and 
the HWQM of the map (outer white line). The three horizontal 
lines indicate the position of cuts across the jet shown in Figure [6] 
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Fig. 6. — Comparison of the normalized synchrotron intensity 
across the jet beam. The overall best model A is indicated by 
heavy thick lines. For comparison model C, which has the most 
pronounced limb-brightening, is shown with thin lines. We plot- 
ted cuts across the jet beam at three different positions along the 
jet, i.e. close to the core (solid/black line), at intermediate dis- 
tance (dashed/red line), and at large distance (dotted/blue line). 
See Figure [5] for the location of the cuts. The lateral coordinate 
is normalized to the width of the jet R a as defined by the sepa- 
rating fieldline ^ a and the intensity is normalized to unity at the 
maximum. 



3.1. Calculation of synchrotron maps 
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Fig. 7. — Comparison of the opening angle calculated from syn- 
thetic maps and the observational data for M87. The lines show 
opening angle profiles as given by the HWQM contour -R ma p of 
the overall best model A (solid line), and model B (dashed line) 
which fits the opening angle data best. Models C (dotted) and 
D (dot-dashed) are shown for comparison. Various symbols repre- 
sent observational measurements. The data points marked by filled 
circles were taken into account in the fitting procedure. The inner- 
most and outermost measurement (open circles) were disregarded 
as explained in the text. 

The calculation of the radio emissivity is done ac- 
cording t he relat i vistic gene ralization of expres sions pre- 
sented by iLaind (|1981f ) and iPacholczvkl (|1970f ) . It is as- 
sumed that the radiating region is optically thin with a 
uniform and isotropic distribution of electrons 

N'(E') oc £'-( 2 «<=+i) (1) 

giving rise to radiation with a spectral index a c —l, i.e., 
S' v oc i/ -Q ° . Note, that primed quantities are measured 
in the comoving frame of the plasma, while unprimed 
quantities refer to the lab frame or are independent of 
the frame of reference. 

We assume further, that the fraction of electron 
number density to proton number density is constant 
throughout the emitting volume. Then the number den- 
sity of electrons n e — J N(E)dE is proportional to the 
density of the plasma, i.e. n e oc p in every frame of 
reference. 

If the radiating plasma moves at relativistic speeds v, 
i.e. T = 1/(1 — (3 2 ) 1 / 2 1 with (3 — v/c, it is convenient 
to evaluate the synchrotron emission in the comoving 
frame X', instead of the lab- or observer frame E. Note 
however, that the velocity of the comoving frame is not 
constant within the emitting region, unless the velocity 
field v is homogeneous, which is not true in our model. 
The synchrotron emissivity in the comoving frame e' is 
related to the magnitude of the magnetic field perpen- 
dicular to the line-of-sight and the electron density as 

e' <xp'\B' xn( os r+\ (2) 

where p' ,B' and n[ os are the electron density, magnetic 
field and line-of-sight unit vector in the comoving. These 
are given in terms of the observers frame quantities 
through the Lorentz-transformations as, 

P = P/T, (3) 



B ' = f B + fT~i^ {v - Bl E ' = (4) 

h{ os =Vn los -(V + l)-^—-. (5) 
1 + 1 c 

We have introduced the Doppler factor T> — (T(l — v- 
nios/c))^ 1 for a compact notation and exploited the fact, 
that in ideal MHD Ohm's law holds as E = —v x B in 
both frames of reference. 

The emissivity in the lab frame e appears Doppler- 
boosted as 

e = V a ° +2 e'. (6) 

The amount of relativistic beaming strongly depends on 
the line-of-sight angle, i.e. the angle between the jet 
axis and the direction to the observer, which we fix to 
#ios = 40°. Finally, the flux in the plane of the sky, /, is 
given by integration along the line-of-sight 

I = J ed£ los . (7) 

The synthetic synchrotron maps are convolved with 
a Gaussian beam to qualitatively match the finite reso- 
lution of observed maps. However, since our synthetic 
maps span more than four orders of magnitude in dis- 
tance from the core, the width of the Gaussian beam 
is not kept constant. Typical radiomaps have a spatial 
resolution corresponding to a couple of observing beams 
across the width of the jet. We therefore use at each 
distance along the jet a Gaussian convolution kernel of 
1-sigma width equal to a tenth of the jet radius at that 
distance, i.e. a(Z) = R a (Z)/ 10. 

The measu r emen ts of the opening angle collected by 
iBiretta et al.l (|2002l ) typically define the jet radius at dis- 
tance Z as the half-width at quarter-maximum (HWQM) 
for the emission across the jet. We adopt this definition 
and refer to it as the jet width of the radiomap 2i? map , 
with 

7(i? map , Z) = max(I(R, Z))/4. (8) 

3.2. The best model 

The synthetic synchrotron emission maps shall quali- 
tatively reproduce two observational constraints: (i) the 
opening angle defined through HWQM (half-width at 
quarter-maximum) in terms of small % 2 values, and (ii) 
the pronounced limb-brightening in terms of small val- 
ues for the mean intensity on the axis over intensity on 
the limb, i.e. (/axis/^iimb)- As a secondary criterion, we 
favor models showing some degree of enhanced emission 
close to the nominal position of HST-1 at 70 pc. 

We calculated synthetic synchrotron emission maps 
and evaluated them according to our two criteria. For 
each of the two criteria we assigned a score between 0..1 
from a sorted list of values for % 2 and (/axis/^iimb), re- 
spectively, and added those to obtain the total score. A 
few formally high scoring models were discarded because 
they did not show clear sign of increased brightness at 
distances 70-200 pc, that could be identified with the 
knot HST-1. The first model satisfying all criteria will 
be referred to as the best model or simply model A. 

The jet in the best model A is launched with a bulk 
velocity of Tj = 2.74 and temperature T = 2.85 mc 2 . 
The separating field line threads the launching surface 
radially at an angle a = 16° with an angular velocity 
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ujj = 3 r g /c. The outer cold disk- wind is launched with 
velocity T d = 1.02. 

For the best model A we plotted spatial 2-dimensional 
synthetic synchrotron emission maps (Figure ©, the in- 
tensity measured on the jet axis as a function of distance 
from the core (Figured]), and the intensity profiles across 
the jet at various positions along the jet axis (Figure 0. 

In Figure [4] we plot the trend along the axis for the 
best model A and compare it with model C, which shows 
the highest local brightness enhancement at the position 
of HST-1 in our sample. The intensity along the axis 
is well described by a power-law in projected distance. 
Both models have similar power-law indices ~ —2.5. 
The model jets dim out faster than the observed jet 
whose power-law index can be estimated to ~ —1.5 from 
the contrast of several published intensity maps. How- 
ever, this discrepancy is not surprising as our model does 
not take into account any micro-physics in order to re- 
energize the electron distribution. Both models show a 
rise of luminosity around 70 pc over the local power-law 
by a factor of ~ 3 and ~ 10, respectively. The location 
of this bright spot coincides with the location of strong 
recollimation shocks, where density and magnetic field 
strength increase. 

Figure [5] shows the synthetic synchrotron map. In or- 
der to increase the contrast we divided the map by the 
trend along the jet axis, i.e. I(R, Z)/I(0, Z). The jet 
beam is well defined, showing large opening angles close 
to the core and becoming almost conical at large dis- 
tance. At the position of the recollimation shock the jet 
width decreases and forms a visible neck. 

The synthetic synchrotron maps show strong gradients 
of intensity across the jet. These are more clearly visible 
in Figure El were cuts across the jet are shown. For 
comparison we show also profiles for model C which has 
the most pronounced limb-brightening. 

The best model A shows limb-brightening already at a 
distance ~ 0.2 pc from the core. In the limb-bright region 
(-faxisMimb) = 0.76 on average before convolution. After 
convolution, the ratio rises to typically 0.85. Model C bi- 
furcates at ~ 0.5 pc and has (/axisMimb) = 0.67 and 0.80 
before and after convolution, respective ly. These values 
agree with those from observation s by [Lv et al.| ((200' 
~ 0.63. iKovalevI (|2008h ~ 0.6, and lKovalev et alJ pOO 
and private communication) 0.6 — 0.8, depending on jet 
region and resolution. We note however, that the limb- 
brightness, as well as the bifurcation distance depend on 
the details of convolution as discussed later on. 

Figure [7] compares the opening angle derived from the 
synthetic maps R map with the observations. Model B 
fits the observations best with % 2 = 5.7, but fails to 
reproduce the limb-brightening. The best model A has 
X 2 = 14.6, which we consider acceptable, in particular, 
if one considers the figure. 

Finally, we compare the width of the jet from the MHD 
models, i? Q , with the jet width of the synthetic syn- 
chrotron maps, i? map - Both curves are superimposed on 
the synthetic map in Figure [5] For many models in our 
sample these two curves are virtually indistinguishable. 
However, for our overall best model, these two curves 
are noticeably different, even if they run almost parallel. 
For this particular model R a underestimates the open- 
ing angle by ~ 20% for all the length of the jet, i.e. 
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Fig. 8. — Synchrotron emissivity in the mid plane of the jet. To 
increase contrast we normalized each row at distance Z by the max- 
imum emissivity across the jet at that particular distance. Close to 
the core the whole jet body contributes to the total emission along 
the line-of-sight. Further down the flow, most of the synchrotron 
emission originates from a thin shell at the outer edge of the jet, 
causing the bifurcation of the jet from center-bright to limb-bright. 
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1.2R a . For other models in our sample, the 



difference is typically not more than 30% — 40%. 

3.3. The origin of limb-brightening 

The synthetic emission map of the best model A in 
Figure [5] shows clear limb-brightening from a distance of 
~ 0.5 pc up to the location of the recollimation shock at 
~ 70 pc and to a lesser extend even beyond. This can 
be seen more clearly in Figure which shows the syn- 
chrotron emissivity e' in the mid plane of the jet. It is 
worth noting, that close to the origin the whole jet body 
emits synchrotron radiation. Further down the jet - and 
certainly beyond 0.2 pc projected distance from the core 
- the emission is dominated by a thin shell at the outer 
edge of the jet, while the inner region close to the axis 
remains relative dark by a factor of more than a hun- 
dred. Within our model this is easily explained noting 
that the synchrotron emissivity equation @ is roughly 
proportional to p' B' 2 . Fieldlines within the relativistic 
jet start out radially from the launching surface, but are 
soon deflected by the outer cold disk wind. This leads 
to strong concentration of magnetic flux, and therefore 
high field strength in a narrow region at the interface 
between both components of the MHD model. The si- 
multaneously occurring increase of density plays only a 
minor role. 

Further down the flow, at around 60 pc, the emissiv- 
ity within the jet body increases again significantly and 
edge-brightening is less pronounced. This region coin- 
cides well with the location of the recollimation shock, 
where the topology of the magnetic field and the distri- 
bution of plasma across the jet change significantly again. 
Beyond the recollimation shock, the emissivity is almost 
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homogeneous across the jet. 

However, the full picture is more complicated 
than that, since the orientation of the highly non- 
homogeneous velocity field relative to the line-of-sight 
is at least equally important. Not only does the poloidal 
velocity vary from one magnetic flux surface to the next, 
but the jet plasma is also rotating. Together this makes 
the Doppler factor and the angle to the line-of-sight in 
the comoving frame highly variable even along the cir- 
cumsphere of the emitting plasma shell. Cross sections 
perpendicular to the apparent axis of the jet reveal, that 
the emissivity along the line-of-sight may be highly asym- 
metric with respect to fore- and background halves of the 
jet. 

4. DISCUSSION 

Limb-brightening strongly constrains the parameters 
of MHD models for the jet of M87. It is not only diffi- 
cult to find parameters resulting in a pronounced limb- 
brightening in terms of low (iaxis/^iimb), but also repro- 
ducing the spatial extend of the limb-bright region. De- 
pending on frequency and resolution, radiomaps for M87 
show limb-brightening from large distances almost right 
down to the core ([Cheung et al.1 120071 : iKovalev et al.1 
2007). However, the bifurcation distance, i.e. the lo- 
cation where the jet morphology transits from center- 
bright to limb-bright, seems to depend strongly on the 
frequency of a particular observation. It is not clear 
a priory to what extent this is due to intrinsic differ- 
ent emission properties at different observing frequen- 
cies, or due to spatial resolution effects. Y. Y. Ko- 
valev (private communication) has performed an anal- 
ysis of t he deep 15 GHz VL BA image of the inner jet 
in M87 (jKovalev et al.ll2007f ). The distance from the 
core at which the bifurcation becomes detectable changes 
from 5 mas at the original 15 GHz resolution to more 
than 100 mas for a tenfold larger beam corresponding 
to 1.5 GHz. We conclude, that variable resolution can 
in principle account for different bifurcation distances at 
different frequencies. 

Limb-brightening is often explained by stratified jets 
consisting of a f ast spine and a slow sheath (see e.g. 
lAlov et al.|[2000l: iGhisellini et al.l 120051 ). In those mod- 
els the slow sheath brightens up in comparison to the 
spine because of Doppler boosting into the 1/L cone. 
For the given large line-of-sight angle of the M87 jet, 
the fast spine is seen from well outside of the Doppler- 
boosted cone and its emission deboosted. Here we sug- 
gest an alternative explanation for limb-brightening in 
AGN, in particular M87. In our model the limb bright- 
ens up due to the concentration of magnetic flux at the 
interface between the relativistic jet and the confining 
non-relativistic disk-wind. The synchrotron emissivity 
is therefore already intrinsically higher in the comov- 
ing frame and limb-brightening does not rely on rela- 
tivistic aberration alone. Incidentally, as seen in Fig- 
ure [3] in our models the spine may at times be mov- 
ing slower than plasma further out which will counteract 
limb-brightening for large inclination angles. 

The exact value of the line-of-sight angle for the jet of 
M87 is still under debate. In this study we have assumed 
an angle 9\ os = 40° which brings it into the upper end of 
the suggested range of values. In order to study the ef- 
fect of the line-of-sight angle on limb-brightening we have 



calculated synthetic maps for model A under an angle 
#ios = 30°. In comparison to the original map, the bifur- 
cation distance moves out to 0.5 pc, limb-brightening is 
more pronounced, and the recollimation shock is brighter 
with regard to the trend along the axis. At the same 
time this obviously enlarges the apparent opening angle 
along the jet and brings the recollimation shock closer 
to the core in disagreement with observations. The sim- 
plest way to bring the opening angle curve back down 
and move the recollimation shock out, is to decrease the 
model's initial opening angle a. However, models with 
lower a tend to show weaker limb-brightening and dim- 
mer recollimation shocks (if present at all). It may there- 
fore be challenging to construct models that reproduce 
all criteria for smaller inclination angles. 

Under the assumption that jet and counterjet are in- 
trinsically identical, one can in principle constrain the in- 
clination angle by measuring the jet-to-counterjet bright- 
ness ratio I + 1 /_ given by 



1 + COS 6>i os 
1 COS 6»i os 



ao+2 



(9) 



where f3 = v/c is the velocity in units of the speed of 
light and a e the ele ctron distribution power-law index. 
IKovalev et al.l (I2007D estimated this ratio to be of the or- 
der 10 - 15, while lLv et all (|2007h measured I+/I- = 14.4 
and constrained the bulk flow velocity to (3 ~ 0.6 — 0.7. 
We calculated a synthetic synchrotron map for the best 
model A pointing in the opposite direction, i.e. the coun- 
terjet at 6>i os = 220°, and estimated I+/I- ~ 13 from the 
on-axis intensities close to the core on either side of it (see 
Fig [9]). The agreement is remarkable, in particular tak- 
ing into account, that one would naively expect a much 
higher value (> 200) from the model's initial Lorentz fac- 
tor Ij = 2, 74, i.e. 0j = 0.93, given in table [1] However, 
eq. (J9|) assumes a homogeneous velocity field along the 
axis and identical intrinsic emissivities on either side of 
the core. Stratification of the jet beam and rotation of 
the plasma about the axis, as is the case in our mod- 
els, can therefore lead to serious misinterpretation if not 
taken into account. 

While fitting the openin g angle to th e obse rvational 
measurements collected bv iBiretta et al.1 (|2002j ). we dis- 
regarded the outermost data point close to the optical 
knot A at 900 pc. All our models fail to reproduce this 
measurement and disregarding it made the least-squares- 
fit more reliable and robust for the remaining data points. 
However, we argue, that at such large distances the in- 
teraction between the environment and the jet, which 
cannot be taken into account by our numerical MHD 
solver, begins to dominate the jet dynamics and thus 
renders our model insufficient. We have also disregarded 
the innermost data point of 60°. This old m easurement 
has been challenged recently. iLv et al.l (|2007h state that 
the opening angle at the jet base is large r than 15 ° 
and therefore consistent with IBiretta et alj ()2002f ) old 
measurement, but than interestingly continue to repeat 
arguments t hat it might be an over estim ation. Maps 
presen ted bv iKrichbaum et all (|2008l ) and IWalker et al.1 
(2008) seem to leave room for smaller opening angles at 
the base. In general, the initial opening angle measured 
from our synthetic emission maps can easily be higher 
than 40°. 
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Fig. 9. — Synthetic synchrotron map for the inner jet of M87 including the counterjet convolved with a constant Gaussian beam of 
cr = 0.03pc. Unlike Fig0this map has not been divided by the trend along the axis, but shows the full dynamic intensity contrast of 50000 
on a logarithmic scale. The contour lines are chosen as indicated on the colorbar in order to highlight the limb-brightening. The counterjet 
is under-luminous by factor of ~ 13. This map is meant to be compared to Fig 1 in Kovalcv ct al. (2007). 



Most of our models that reproduce the measured open- 
ing angle as a function of distance as well as the limb- 
brightening, do at the same time exhibit a recollimation 
shock 1 showing some degree of locally enhanced emission. 
We identify this recollimation shock with the bright knot 
HST-1. In most of our models the emission at the bright 
knot does not increase by more than a factor of 10 com- 
pared to the general trend, while for the real HST-1 the 
increase in emission level is probably higher. However, 
this is not surprising, since our MHD model does not in- 
clude any micro-physics as e.g. particle acceleration, that 
could be responsible for further enhancing the emission 
level. It is still noteworthy, that the conditions in the 
vicinity of the recollimation shock of the best model A, 
are similar to those postulated from SSC m odels in order 
to fit t he TeV emission in M87. In particular lHarris et al.l 
(2003) recently presented a SSC model assuming, that 
the TeV emission was originating from HST-1. In our 
model the Doppler factor in the vicinity of the recollima- 
tion sh ock is rather l ow T> ~ 1 — 2. For those Doppler 
factors, lHarris et alj (|2003[ ) calculated a magnetic field 
strength of B' = (13 — 3.7) mG, which matches very well 
B' ~ 10 mG in our model. 

5. SUMMARY AND CONCLUSIONS 

We have calculated self-consistent global MHD models 
and synthetic optically thin synchrotron maps with the 
aim to reproduce the reported opening angle distribution 
as well as the morphological structure, in particular limb- 
brightening, of the jet of M87. We have applied different 
criteria to quantify the agreement between models and 
observations. All criteria can be satisfied to high degree 
by individual models in our database. However, no model 
does satisfy simultaneously all criteria exceedingly well. 
We identified a best model as compromise between the 
opening angle distribution and the high degree of limb- 
brightening. This model satisfactorily meets our criteria, 
b ut does not p e rform best in any single one. 

iGracia et al.l (|2005f ) used the shape of the fieldline sep- 
arating the two regions on the boundary surface to define 
the opening angle. As shown in this work, this particular 
fieldline and flux contours of synthetic maps are paral- 
lel for a wide range of models and distances along the 
jet. However, the first method may under-estimate the 
opening angle for some parameters by typically 30%. 

The morphological structure across the jet can be re- 
produced in principle. A wide range of models show limb- 

1 We also have a small range of models, that reproduces the 
opening angle wit hout showing a strong recollimation shock (see 
IGracia et al. 2005], Figure 2), but those models do in general not 
show clear limb-brightcning. 



brightening away from the core. The transition from 
center-bright to limb-bright, the bifurcation of the jet, 
generally occurs well below 1 pc projected distance. A 
quantitative comparison with specific observations, how- 
ever, requires careful modeling of the convolving beam 
and may require including opacity effects, i.e. optically 
thick emission. Both effects are frequency dependent. 

In our models limb-brightening is due to a strong in- 
crease of the comoving frame synchrotron emissivity in 
a relatively thin shell near the outer edge of the visible 
jet as result of a concentration of magnetic flux at the 
interface between the relativistic outflow and the non- 
relativistic disk-wind, i.e. the environment. In contrast 
to other models limb-brightening here is not an immedi- 
ate result of deboosted emission in a slow moving sheath; 
in fact the 'sheath' in our models typically moves slightly 
faster than the jet body. 

The optical knot HST-1 has received much attention. 
Most MHD models with acceptable fits to the opening 
angle distribution feature a shock at the distance 70-200 
pc. This is particularly true for models with large ini- 
tial opening angle a. The shock is due to the fact, that 
the jet is not in lateral force equilibrium at the launch- 
ing surface and over-expands. Later, it is forced back 
towards the axis by magnetic hoop stress and prod uces 
a recollimation shock (jBogovalov fc Tsinganosll2"005[ ) . In 
all cases the shock is visible as a bright spot in the syn- 
thetic maps. 

We performed our parametric studies unbiased towards 
the existence of a recollimation shock. All high ranking 
models have a recollimation shock. Only after ranking 
our models, we favored for illustrative purposes models 
with pronounced brightness increase at some point along 
the axis. We therefore conclude, that the optical bright 
knot HST-1 is a general feature of all models matching 
our selection criteria and cannot be seen independently 
from them. In particular, we suggest, that any physical 
model explaining the opening angle and morphological 
structure will simultaneously also account for HST-1. 

In addition our best fitting model is consistent with 
a number of observational constraints such as the mag- 
netic field in the knot HST-1, and the jet-to-counterjet 
brightness ratio. 
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